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The site percolation threshold for the random Voronoi network is determined numerically for 
the first time, with the result p c = 0.71410 ± 0.00002, using Monte-Carlo simulation on periodic 
systems of up to 40000 sites. The result is very close to the recent theoretical estimate p c m 0.7151 
of Neher, Mecke, and Wagner. For the bond threshold on the Voronoi network, we find p c — 
0.666931±0. 000005, implying that for its dual, the Delaunay triangulation, p c = 0.333069±0.000005. 
These results rule out the conjecture by Hsu and Huang that the bond thresholds are 2/3 and 1/3 
respectively, but support the conjecture of Wierman that for fully triangulated lattices other than 
the regular triangular lattice, the bond threshold is less than 2sin7r/18 s=s 0.3473. 

PACS numbers: 



I. INTRODUCTION 

The Voronoi diagram [T] for a given set of points on a 
plane (Fig.JTJ is simple to define. Given some set of points 
P on a plane K 2 , the Voronoi diagram divides the plane 
R 2 into polygons, each containing exactly one member of 
P. Each point's polygon cordons off the portion of K 2 
that is closer to that point than to any other member of 
P. More precisely, the Voronoi polygon around Pi 6 P 
contains all locations on K 2 that are closer to than to 
any other element of P. The total Voronoi diagram is the 
set of all the Voronoi polygons for P on R 2 ; the Voronoi 
network is the set of vertices and edges of the Voronoi 
diagram. 

The dual to the Voronoi diagram is interesting in its 
own right. Known as the Delaunay triangulation |2J (see 
Fig. it can be defined independently of the Voronoi 
diagram for the same set of points P on I 2 : it is simply 
the set of all possible triangles formed from triples chosen 
out of P whose circumscribed circles do not contain any 
other members of P (Fig. [3]). The Delaunay triangulation 
and the Voronoi diagram for the same set of points can 
be seen in Fig. [4] Note that while the members of P are 
sites in the Delaunay triangulation, they are not sites 
in the Voronoi network, whose sites are the vertices of 
the polygons; also note that the edges of the Voronoi 
diagram lie along the perpendicular bisectors of the edges 
of the Delaunay triangulation — however, the edges of 
the Voronoi diagram do not always intersect the edges 
of the Delaunay triangulation, as seen in Fig. [4] The 
Delaunay triangulation represents the connectivity of the 
Voronoi tessellation of the surface. 

There are many algorithms for constructing these net- 
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FIG. 1: Voronoi diagram with a Poisson distribution of gen- 
erating points. 



works. The fastest ones run in 0(n log n) time for general 
distributions of points HJ EJ E], where n is the num- 
ber of generating sites, and this has been proven to be 
the optimal worst-case performance 0]. For a Poisson 
distribution of points on the plane, there are many O(n) 
expected-time algorithms [7] |9j [10] . 

In addition to being theoretically interesting [TT1 fT51 
[T3j [14l [T5] , both the Voronoi diagram and the Delaunay 
triangulation are widely used in modeling and analyz- 
ing physical systems. They have seen use in lattice field 
theory and gauge theories [16] . analyzing molecular dy- 
namics of glassy liquids [17], detecting galaxy clusters 
|18j . modeling the atomic structure and folding of pro- 



FIG. 2: Delaunay triangulation for the same set of generating 
points as in Fig.[TJ The generating points become the vertices 
in this network. 

teins [THl [20] , modeling plant ecosystems and plant epi- 
demiology |21j . solving wireless signal routing problems 
[221 E3] , assisting with peer-to-peer (P2P) network con- 
struction [21], the finite-element method of solving dif- 
ferential equations [25] , game theory [25] [27J , modeling 
fragmentation [21], and numerous other areas [2U1 I3U] . 

Percolation theory is used to describe a wide variety of 
natural phenomena [3JJ [32] • In the nearly seventy years 
since the first papers on percolation appeared [33] [34], 
it has become a paradigmatic example of a continuous 
phase transition. For a given network, finding the crit- 
ical probability, p c , at which the percolation transition 
occurs is a problem of particular interest. p c has been 
found analytically for certain 2D networks [35] [351 I3T] ; 
however, most networks remain analytically intractable. 
Numerical methods have been used to find p c for many 
such networks, e.g., [3TJ [38] [39] HQ] SU S2J H3] . 

In this paper we consider the percolation thresholds 
of the Voronoi and Delaunay networks for a Poisson dis- 
tribution of generating points, as represented in Figs. [TJ 
[2] and [4] There are four percolation thresholds related 
to these two networks: the site and bond percolation 
on each. Being a fully triangulated network, the site 
percolation threshold of the Delaunay network is exactly 
p stie,Dei _ i [33J[2U[25J[2B]. This result has recently been 
proven rigorously by Bollobas and Riordan [UJ. Some- 
what surprisingly, a search of the literature revealed no 
prior calculation of the site percolation threshold of the 
Voronoi network at all, despite the widespread use of 
such networks. A prediction for its value has recently 
been made by Neher, Mecke and Wagner [48]; they use 
an empirical formula to predict p<^ te > Vor — 0.7151, but 
they too were unable to find any previous calculation 



FIG. 3: The Delaunay triangulation for a set of five points, 
along with the associated circumcircles. 




FIG. 4: The Delaunay triangulation (dotted lines) superposed 
on the Voronoi diagram (solid lines), its dual graph, for a set 
of Poisson-distributed generating points. 



of this value, either analytically or numerically. (There 
are places in the literature (e.g. [29] ) where the Voronoi 
"site" threshold is listed as 1/2. This is true if the "sites" 
are taken to be the generating points from which the dia- 
gram is created, rather than the vertices of the diagram. 
Thus, this is actually the Voronoi tiling threshold, i.e. 
the percolation threshold of the Voronoi polygons, which 
is in turn equivalent to the Delaunay site threshold, well 
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known to be 1/2.) 

The bond thresholds of the Voronoi and Delaunay net- 
works are complementary, 

bond 7 Vor i bond, Del /i\ 

Pc 1 Pc ' V / 

because these networks are dual to one another [15] . The 
first numerical measurement of the bond threshold for 
either network seems to be that of Jerauld et al. [ST)] , 
who in 1984 found p^Dei = g.332. Shortly there- 
after, Yuge and Hori |51j performed a renormalization 
group calculation which yielded pbond.Dei _ q.3229. In 
1999, Hsu and Huang [SJ] found p b c ond - D ^ = 0.3333(1) 
and p b c ° nd ,Vor = 0.6670(1) through Monte Carlo meth- 
ods. (The numbers in parentheses represent the errors 
in the last digits.) These values led them to make the 
intriguing conjecture that the thresholds are exactly 1/3 
and 2/3 respectively. There is, however, no known the- 
oretical reason to believe that this conjecture is true. In 
order to test this conjecture, and to find the site perco- 
lation threshold of the Voronoi network, we have carried 
out extensive numerical simulations, as detailed below. 
In section [TTJ we describe our methods, and in section |III| 
we discuss our results and compare them to the thresh- 
olds of several related lattices, and also discuss the cov- 
ering graph and further generalizations of the Voronoi 
system. Conclusions are given in section |TV| 

II. GENERATING ALGORITHMS AND 
ANALYSIS TECHNIQUES 

A. Delaunay /Voronoi generation algorithm 

In order to avoid edge effects in the networks when 
growing percolation clusters, and to make it possible to 
use m ore sites as seeds for those clusters (see subsection 
IIB| , we wished to create Voronoi and Delaunay net- 
works with periodic boundary conditions. While popular 
fast algorithms for generating Voronoi and Delaunay net- 
works exist, most notably the Quickhull algorithm [53] . 
these do not generally support periodic boundary con- 
ditions. We therefore created our own fairly straightfor- 
ward algorithm for generating the desired networks. Af- 
ter coming up with it independently, we later found that 
it falls into the class of expected-linear-time algorithms 
known as Incremental Search [8]. The basic outline of 
the algorithm is as follows: 

1. Divide the region in which the generating points 
(vertices of the Delaunay triangulation) are located 
into squares of equal size ( "bins" ) . 

2. Find a single Delaunay edge by picking a point at 
random and searching through its bin and neigh- 
boring bins to find the point which is its nearest 
neighbor. 

3. Given a Delaunay edge and a "side" to look on (im- 
mediately above or below the edge), determine the 



third point in the Delaunay triangle by looking at 
the radii of the circumcircles of the triangles formed 
by that edge with each point in its bin and all of 
the neighboring bins. 

4. Look at the other Delaunay triangles that have 
been found in that bin and neighboring bins to 
make sure this new triangle is not a duplicate of 
one that has already been found. If it is not, find 
which of its neighbors have already been discovered 
and mark them as its neighbors, and vice versa. 

5. From the list of neighbors of the current triangle, 
figure out which of its edges are not already shared 
with neighbors (if any), and, if there are any un- 
shared edges, whether the missing neighbor should 
be above or below the edge. 

6. Repeat steps 3-5 until there are no unprocessed 
edges left, at which point the Delaunay triangu- 
lation is finished. Because the neighbors of each 
triangle are known, this algorithm also yields an 
adjacency list of the sites on the Voronoi diagram 
(because the Voronoi diagram is dual to the Delau- 
nay triangulation). 

This algorithm is significantly easier to implement with 
periodic boundary conditions, because every triangle is 
guaranteed to have exactly three neighbors. Further- 
more, the imposition of periodic boundary conditions 
also gives the Delaunay and Voronoi networks a very use- 
ful property: there are always exactly twice as many De- 
launay triangles (Voronoi sites) as there are generating 
points (vertices of the Delaunay network or polygons in 
the Voronoi network) for a given diagram. This is a con- 
sequence of the more general fact that that the number of 
faces (triangles) must be double the number of vertices 
(sites) for any fully triangulated network with doubly- 
periodic boundary conditions in two dimensions. This 
fact follows from Euler's formula and is proven in the 
appendix. This simple relation makes it easier to spot 
certain kinds of errors in the code, because improperly 
written code is rather unlikely to consistently produce the 
proper number of sites for the given number of generat- 
ing points. Using this algorithm, we generated thousands 
of Voronoi diagrams of 40,000 sites each. Fig. [5] shows 
an example of a smaller Delaunay triangulation created 
with this algorithm. 



B. Percolation cluster growth and finding p c 

The Leath-type epidemic growth method [54] that we 
used involves growing a large number of percolation clus- 
ters in order to find p c . For site percolation clusters, we 
start with a seed site somewhere on the network. Each of 
its neighbors is turned on with probability p or off with 
probability 1 — p. Neighbors of active sites are then vis- 
ited and the procedure is repeated for all their previously 
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FIG. 5: (Color online) A Delaunay triangulation with periodic 
boundary conditions (i.e., on a torus), created from n — 300 
generating points. Because the surface has periodic boundary 
conditions, there are exactly 2n — 600 triangles here. Note 
the corresponding shapes of the outline on opposite edges, 
because this diagram has been "unrolled" from a torus. 



unvisited neighbors; the cluster either dies out naturally 
or is stopped by the program when it hits a cutoff size of 
1000 sites. For bond percolation clusters, an analogous 
algorithm is used in which sites are simply never turned 
off, as it is the bonds between sites that are pertinent. 

Due to the fact that our Voronoi diagrams are finite 
and are generated from a random Poisson distribution of 
points, each diagram yields a slightly different effective 
value of p c ; therefore, we had to generate many diagrams. 
This could have been very computationally expensive, 
but the choice of periodic boundary conditions helped 
here as well. Because there are no edges to the diagrams, 
we were able to place the seed point for a cluster at any 
site on a diagram, rather than being limited to a small 
subset of sites near the center. This meant we were able 
to use many widely separated seed points to grow clus- 
ters on each diagram, which reduced the impact of each 
seed point's immediate neighborhood upon the value of 
p c obtained for each diagram. This, in turn, dramati- 
cally reduced the number of distinct diagrams we needed 
to obtain a particular level of precision. Specifically, we 
grew 8 x 10 5 clusters of up to 1000 sites on each of 800 
diagrams, for a total of 6.4 x 10 s clusters grown at each 
value of p. We then repeated this process at each of 
various p near p c to generate the plots in the next sec- 
tion. Finally, this process was done twice — once for site 
percolation and once for bond percolation, both on the 



Voronoi network. 

Because the percolation clusters are cut off before they 
can become large enough to wrap around the network, 
the clusters effectively see the diagram as infinite in size. 
Thus, their size distribution can be used to obtain an un- 
biased estimate of P s , the probability that a percolation 
cluster will grow to be at least size s (for s < 1000) on an 
infinite network. At the critical threshold p c , P s ~ s 2 ~ r 
as s — > oo, where r = 187/91 for the two-dimensional 
percolation cluster universality class [31]. (It is expected 
that the critical exponents here are the same as for regu- 
lar two-dimensional lattices.) In the scaling region, where 
s is large and p — p c is small such that s a (p — p c ) is con- 
stant (with a = 36/91), P s behaves as 

P s ~ As 2 -Tf(B(p- Pc )s°), (2) 

where A and B are non-universal metric constants spe- 
cific to the system being considered, and f{x) is a univer- 
sal scaling function analytic about x = 0. If we operate 
close to p c such that B(p — p c )s a -C 1, then we can make 
a Taylor-series expansion of f(x) to find 

P 8 ~s 2 - T (A + D(p- Pc )s a + ...), (3) 

where D is another constant. Thus, plotting C s = 
P s s T ~ 2 vs. s a should yield a straight line at large s 
when p is near p ci and that line will have a slope of zero 
when p — p c . Fig. [6] shows several such plots for site 
percolation clusters on the Voronoi network, and Fig. [7] 
shows several plots for bond percolation clusters on the 
same. C s does indeed approach a linear function for large 
s in these plots, albeit far more quickly for bond perco- 
lation than for site percolation, with p^ tte , Vor ~ 0.7141 
and p bond,Vor _ 0.66693. 

Unfortunately, for smaller s there are deviations in C s 
due to finite-size effects, and these are quite apparent for 
site percolation, even at the largest values of s we were 
able to investigate. Exactly at p c , one expects 

P s ~ s 2 - T (A + Es- n + ...) (4) 

as s — > oo, where £ is a constant and £1 w 0.6 — 0.8 is 
the corrections-to-scaling exponent |55j . Similar devia- 
tions should occur when p is close to p c - In the case of 
site percolation on the Voronoi network, these finite-size 
effects make it difficult to determine when C s has a truly 
horizontal asymptote; thus, it is not possible to use the 
above method to find p^ lte - Vor to much greater precision 
than four digits when s < 1000. 

The most straightforward way to solve this problem 
would be to grow larger site percolation clusters, using a 
larger system to insure that wrap-around does not occur. 
However, because of the computational time that would 
be required to do that, we instead used a more sensitive 
method to find p c that takes the finite-size corrections in 
Q into account. 

Eq. Q implies that, at p c , C s - C s/2 = E(l - 2 n ) S - n 
to leading order. This means it's possible to estimate 
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FIG. 6: (Color online) Epidemic site percolation cluster 
growth on the Voronoi diagram, for p = 0.71407, 0.71409, 
0.71411, and 0.71413, from bottom to top on the right. 
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FIG. 7: (Color online) Epidemic bond percolation cluster 
growth on the Voronoi diagram, for p = 0.66691, 0.66693, 
0.66695, 0.66697, and 0.66699, from bottom to top on the 
right. Note that these plots approach their linear asymptotes 
far more rapidly than those for site percolation clusters, as in 
Fig. [6] also note the difference in the vertical scale between 
the two figures. 



FIG. 8: (Color online) A zoomed-in portion of Fig. [7] The de- 
viations from the horizontal in the asymptotes for each curve 
can be seen more clearly here. (The y-axis values are given 
from a reference of C 3 = 1.058). The least-squares linear fits 
for the curves are also on this plot. 
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FIG. 9: (Color online) {l est at p = 0.71407, 0.71409, 0.71411, 
0.71413 (top to bottom on right) for site percolation on the 
Voronoi diagram. 



III. RESULTS AND COMPARISON WITH 
RELATED LATTICES 



directly from |55j 



nr' = -b 8a ( c ' c ;! 2 )- (5) 



Thus, in the regime where s is small enough that the 
finite-size effects of Q matter, yet large enough that 
higher-order corrections are unimportant, f2f st should 
approach a constant f2 when p = p c . When p ^ p c , there 
will be deviations due to scaling. Plots of f2f st vs. In s 
for several values of p can be seen in Fig. [9] these yield 
the result for p^- te < Vor found in the following section. 



A. p c for site and bond percolation on the Voronoi 
diagram 

Examining Fig. [9j it can be seen that fif st approaches 
a constant for large s for p w 0.71409 — 0.71411, and we 
conclude 

p f e ' Vor = 0.71410 ± 0.00002 , (6) 

where the error bars are meant to indicate one standard 
deviation of error. This plot also gives us a rough value 
of 0.65 for — close to the value of f2 found for the 
Penrose rhomb quasi- lattice |55j . 

We used the method of plotting C s = P s s T ~ 2 vs. s a , 
outlined in the previous section, to find the bond per- 
colation threshold of the Voronoi diagram. Taking the 
results shown in Figs. [7] and [HJ we see immediately that 
pbond.Vor ^ q. 66693. Because finite-size effects were not 



6 



01 

§~ 0.00002 



0.66690 0.66692 0.66694 0.66696 0.66698 0.66700 

Pc 



FIG. 10: Slopes of the lines fitted in Fig. [5] versus the values 
of p used for each line, along with a best fit line. 



significant for bond percolation, we were able to find ex- 
cellent least-squares linear fits to the asymptotic portions 
of the curves in Fig. [8] By plotting the slopes of these 



lines against the values of p used (see Fig. llOj) , we were 
able to solve for the value of p that would yield a slope 
of zero; this should be p c . This technique yielded a more 
accurate estimate: 



P 



bond, Vor 



0.666931 ±0.000005 



(7) 



which by (0 implies p b c ond - Dd = 0.333069 ± 0.000005. 
We considered various contributions to the stated error. 
First of all, it is unclear precisely where the linear regime 
begins in Fig. [8] and this leads to some uncertainty in 
the slopes we measured from the best-fit lines. Statis- 
tical effects of course are a source of error. However, a 
somewhat larger source of uncertainty turned out to be 
the error involved in reusing the same diagram multiple 
times during cluster growth — even with different seed 
points, there is a distinct likelihood that the same part of 
the non-uniform diagram will be sampled. To estimate 
this error, we considered our usual runs of 800,000 sam- 
ples on 10 different diagrams at p — 0.666931 and looked 
at the variation in the curves of C s vs. s a (Fig. 11). In 
contrast, we also looked at 10 runs of 800,000 samples 
each on the same diagram, to gauge the purely statisti- 
cal error. We found the errors in the previous case larger 
than in the latter. Using the measured standard devia- 
tion 1.5 x 10~ 4 and dividing by -\/800 for the 800 runs we 
actually used in our simulations for each value of p, we 
estimate a final error of ±0.000005 in the slopes of C s , 
as indicated in the error bars of FigJTO] Finally, because 
the slope of the fitted line in Fig. |10| (which equals the 
coefficient D in Eq. ^) is nearly 1, we estimate that the 
final error in p c is ±0.000005. Note that the runs for the 
five values of p were each done on 800 different diagrams, 
so there is no systematic error among the least-squares 
fit lines drawn in Fig. [8j Because of this, we believe our 
error bars are conservative. 

The results for the thresholds are summarized in Table 
UJand discussed further below. 
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FIG. 11: (Color online) Comparison of multiple-diagram and 
single-diagram bond percolation cluster growth on finite peri- 
odic Voronoi networks. Each curve represents 8 x 10 s clusters 
grown on 10 different diagrams (upper) or on 10 identical dia- 
grams (lower), at p = 0.666931. The mean value of the slopes 
are —1.8 x 10" 5 (upper) and —1.1 x 10~ 4 (lower), and the 
standard deviations are 1.5 x 10" 4 (upper) and 4.7 x 10~ 5 
(lower) . 



TABLE I: Results for percolation thresholds of Voronoi and 
Delaunay networks. Numbers in parentheses represent errors 
in last digits. 

"network z pf e p b c ond 

Voronoi 3~ 0.71410(2) 0.666931(5) 
Delaunay 6 (avg.) 0.5 (exact) 0.333069(2) 



B. Comparison with thresholds of related lattices 

The Voronoi network has a uniform coordination num- 
ber z equal to 3. In Table [IT] we compare the site 
and bond thresholds of the Voronoi diagram with sev- 
eral other lattices with z = 3, listed in descending or- 
der of threshold values. In the Grunbaum-Shepard nota- 
tion, (3 a3 ,4 a4 , . . . ) describes a lattice with a 3 triangles, 
d4 quadrilaterals, etc., per vertex. The Archimedean 
lattices (3,12 2 ), (4,6,12), and (4,8 2 ) are illustrated in 
[B*rJl IST1 IS"5] . The martini lattice was introduced in [55] 
and can be represented by (3/4) (3, 9 2 ) + (1/4) (9 3 ). 

We also list in Table III] for each lattice the generalized 
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TABLE II: Thresholds of lattices with uniform coordination 
number z = 3, also showing the filling factor / and polygon 
variance ^/n 2 . a Ref. [56], °Ref. [64], c Ref. [55], d Ref. g5J, 
e this work, / Ref. [40], 9 Ref. [66], *exact. 



lattice 



IP 2 ) 
martini 
(4,6,12) 
(4,8 2 ) 
Voronoi 
honeycomb 



l^ 2 t 



{>,: 



bond 



0.5 0.39067 

0.25 0.47493 

0.222222 0.48601 

0.111111 0.53901 

0.049468 0.57351 

0.0 0.60460 



0.807901*" 

0.764826* c 

0.747806 a 

0.729724 11 

0.71410 e 

0.697040 aJ 



0.740422° 

0.707107* d 

0.693734 c 

0.676802 c 

0.666931 e 

0.652704* 9 



filling factor /, denned as 



E7T 
a n cot — 
n 



n>3 



(8) 



which generalizes Scher and Zallen's definition of / for 
lattices not necessarily composed of regular polygons [BD] . 
The / has been shown to provide a good correlation to 
site percolation thresholds for a variety of lattices. To cal- 
culate / for the Voronoi network, we use 63 = 0.0112400, 
64 = 0.1068454, etc., from [ST], where b n = 2a n /n is 
the fraction of n-sided polygons in the system, satisfying 
En b n = I and n = J2n nb n — 6 for z = 3. 

In Table [U also list the fluctuations in the number of 
the sides of the polygons for each lattice, 



2 —2 

n z — n 



(9) 



which is equal to the fluctuations in the coordination 
number of the dual lattice. It follows from Euler's for- 
mula that the average number of sides of the polygons, 
n = Y^, n nb n , in any 3-coordinated network is exactly six. 
For the Voronoi diagram, /1 as 1.7808116990 is known ex- 
actly as an integral [BTl |6"2"1 R)3"] . 

In Fig. [12] we plot the thresholds given in Table |TT] 
as a function of /. The thresholds fit well to a linear 
relation, as can be seen in the figure. In general, for bond 
percolation, / is not effective in correlating thresholds, 
which depend strongly upon the coordination number z. 
However, for networks with fixed z — 3, we find that the 
correlation of the bond thresholds with / is quite good. 

We can fit the linear behavior of p c (f) using just data 
from exact results, with no numerical input. For site 
percolation, we use the exactly known thresholds for the 
(3,12 2 ) and martini lattices, while for bond percolation 
we use the martini and honeycomb lattice results, and 
find: 



fc 

bond 

Pc 



-0.5116/ + 1.0078 , 
-0.4195/ + 0.9063 . 



(10) 



These equations imply for the Voronoi diagram (where 
/ = 0.57351), p f^ Vor = 0.7143 and p&°«<*,vbr = .6657, 
which are evidently excellent estimates. Thus, the 



thresholds for the Voronoi diagram are consistent with 
other lattices with respect to the filling factor. A similar 
plot of thresholds versus the fluctuations also shows con- 
sistent behavior between the Voronoi results and those 
for these other lattices. 

Finally, the result for the bond threshold for the 
Voronoi network implies the site threshold for the 
Voronoi covering graph, shown in Fig. |13| The cover- 
ing graph (or line graph) for a given network is defined 
as the graph that connects the centers of the bonds to- 
gether, and converts the bond percolation problem on 
that network to a site problem. Thus, p^ te ^ VorCov — 
p bond,Vor _ 0.666931. The covering graph is a kind of 
randomized kagome diagram, consisting of triangles con- 
nected together. Using similar arguments given in |67j 
for generalized kagome lattices, one can find an estimate 
for the bond threshold of the covering lattice, with the 
prediction p boni,VorCov _ q.53618, as well as an estimate 
for the site-bond threshold for the Voronoi diagram. De- 
tails will be given elsewhere |68J. 



IV. CONCLUSIONS 

We have determined the site percolation threshold for 
the Voronoi network, for the first time and to high pre- 
cision, with the result p^e.Vbr = 0.71410(2). We reiter- 
ate that this is not the well-known threshold (1/2) of the 
polygonal tiles of the Voronoi tessellation, which is equiv- 
alent to site percolation on the Delaunay triangulation, 
but rather the threshold for the 3-coordinated diagram of 
all the Voronoi polygons. Our Monte-Carlo result is very 
close to the prediction 0.7151 of Neher, Mecke and Wag- 
ner |48j . and confirms their empirical procedure based 
upon the Euler characteristic. 



0.85 



a" 0.8 

I 75 
I 



0.7 



o 

X) 



0.55 



0.6 



0.55 



(3,12 2 ) 






^martini 




(4,6,12) 








^"^^Vo ronoi 







0.35 0.4 0.45 0.5 0.55 0.6 
generalized filling factor f 



0.65 



FIG. 12: Thresholds vs. generalized filling factor of Eq. Q 
for site (top) and bond (bottom) percolation for the systems 
of Table [n] The lines show least-squares fits to all of the data 
points. 
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FIG. 13: Covering graph of a Voronoi network. 



We also determined the bond threshold of the Voronoi 
network. Our result p b c ° nd > Del = 0.333069(5) is consistent 
with Jerauld et al.'s result 0.332 [50] and close to Hsu and 
Huang's value 0.3333(1) [55], but runs counter to the lat- 
ter authors' conjecture that this threshold is exactly 1/3. 
It is interesting to note that 1/3 is the value predicted 
by the general (approximate) bond-threshold correlation 
p c w d/[(d — l)z] given by Vyssotsky et al. [HH] for z = 6 
and dimension d — 2. 

We made comparisons of our results with thresholds 
of other lattices with the same coordination number 
(z = 3), and found that the Voronoi thresholds are what 
one would expect based upon correlations with the filling 
factor / for both the site and bond problems. 

Wierman has conjectured [70] that 2sin7r/18 w 0.3473, 
the bond threshold of the regular triangular lattice, is the 
maximum possible bond threshold for any fully triangu- 
lated network, and that no other fully triangulated net- 
work has a bond threshold greater than or equal to that 



value. We indeed find that the bond threshold of the 
fully triangulated Delaunay network is consistent with 
this conjecture. 

For future work, it would also be interesting to look at 
thresholds for other random systems, such as Johnson- 
Mehl tessellations [7T] or the graph formed by the ran- 
dom distribution of lines in a plane |72] . Finding thresh- 
olds in Voronoi systems of higher dimensions is another 
interesting open problem. 

APPENDIX: Proof that F = 2V for any fully 
triangulated network with doubly-periodic boundary 
conditions in two dimensions 

We take advantage of the Euler relation for polyhe- 
dra to prove the desired fact about fully triangulated 
networks. A network on a square surface with doubly- 
periodic boundary conditions is topologically equivalent 
to placing the network on the surface of a torus; this net- 
work, in turn, can be seen as a polyhedron on the surface 
of the torus. Thus, the Euler relation for polyhedra ap- 
plies: 



V-E + F 



Xtor 







where V is the number of vertices on the polyhedron, 
E is the number of edges, F the number of faces, and 
Xtoms the Euler characteristic for the 2-torus, which is 
zero. Because every face has exactly three edges (i.e., the 
network is fully triangulated), and every edge is shared 
by exactly two faces (the network has no boundary) , we 
have E — 3F/2, and we can rewrite the Euler relation as 
follows: 



3F 
~2~ 



and thus F 



V 

2V. QED. 



F = V 



F 
2~ 
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